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^— V , We illustrate how the midpoint and iterative cone (with progressive removal) algorithms 

fv^ ' fail to satisfy the fundamental requirements of infrared and coUinear safety, causing 

divergences in the perturbative expansion. We introduce SISCone and the anti-fct 
algorithms as respective replacements that do not have those failures without any cost 
at the experimental level. 
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^^ , The general picture. Jets are an important toolin hadronic physics and theywill play 

^^ ■ a predominant role at the LHC. By defining jets one aims at accessing, from the final-state 

particles, the underlying hard parton-lcvel processes. Since a parton is not a well-defined 
r^ ' object, a jet definition is also not unique. 

Two broad classes of jet definitions exist. The first one works by defining a distance be- 
pi i' tween pairs of particles, performing successive recombinations of the pair of closest particles 

O , and stopping when all resulting objects are too far apart. Algorithms within that class differ 

r^.' by the definition of the distance, frequent choices being dfj = min(A:^j, fc^^)(Ay^^ -I- A^ifj) 

for the kt algorithm [1], and d?- = (Ay? -I- A0? ) for the Cambridge- Aachen algorithm [2]. 
Cone algorithms make up the second class, where jets are defined as dominant directions 
of energy flow. One introduces the concept of stable cone as a circle of fixed radius R in the 
fvj ' y~(t' plane such that the sum of all the momenta of the particles within the cone points in the 

f^ . same direction as the centre of the circle. Cone algorithms attempt to identify all the stable 

^^ ' cones. Most implementations use a seeded approach to do so: starting from one seed for the 

r — . centre of the cone, one iterates until the cone is found stable. The set of seeds can be taken 

f^ ' as the set of initial particles (sometimes over a pt threshold) or as the midpoints between 

00 . previously-found stable cones. As we shall see, this iterative method fails to identify all the 

stable cones, leading to infrared or coUinear unsafety in the perturbative computations. 

Cone algorithms can be split in two subclasses according to how they deal with the 
fact that stable cones may overlap. Cone algorithms with split-merge, identify the hardest 
overlapping pair of stable cones and merge (split) them if they share more (less) than a 
C^ ' fraction / of the hardest cone. JetClu and midpoint are typical representatives of that 

subclass. Cone algorithms with progressive removal start with the hardest unclustered 
particle, iterate from there until a stable cone is found and call it a jet. Its contents are 
removed and one starts again with the remaining particles. The iterative cone is the typical 
example of a subclass with the particular feature that hard jets are fully conical. 

The SNOWMASS accords have established a series of requirements that any jet algorithm 
has to fulfil. These are basically that one can use the algorithm for theoretical computations, 
e.g. it gives finite perturbative results, as well as for experimental purposes, e.g. it runs fast 
enough and has small corrections from hadronisation and the underlying event. 

We illustrate in these proceedings [3] that midpoint and the iterative cone fail to give 
finite perturbative results. We introduce SISCone and the anti-fct algorithms as solutions to 
those problems that do not spoil the experimental usability. 
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Figure 1: Stable cones found by the midpoint algorithm for a 3-particle event (left) and for 
the same event with an additional infinitely soft gluon (right). 



SISCone as a replacement for the midpoint algorithm. Let us consider the 
3-particle event displayed in Fig. [Ua). When clustered with the midpoint algorithm, 2 
stable cones are found, leading to two jets: one with particles 1 and 2 and a second one with 
particle 3. If one adds to that hard event an infinitely soft gluon as shown in Fig. [^b), 
a third stable cone is found and the three hard particles are clustered in a single jet. This 
change in the jet structure upon addition of soft particles, a phenomenon which happens 
with infinite probability in perturbative QCD, gives rise to divergences in the pcrturbative 
expansion and proves that the midpoint algorithm is infrared unsafe. 

This problem arises from the fact that the seeded approach misses stable cones — here 
the one containing particles 2 and 3 in Fig. [Ija). The workaround to restore IR safety 
is thus to find a seedless method that provably identifies all the stable cones. This is 
notoriously complex: a naive approach testing the stability of all subsets of particles [4] has 
a complexity of order 7V2^ for N particles which is much slower than the 0{N^) complexity 
of the midpoint algorithm, making this solution unusable for experimental purposes. 

The solution [6] is to use the geometrical obser- 
vation that any enclosure in the y — (j) plane can be 
moved without changing its contents until it touches 
two points. Browsing all pairs of particles allows thus 
to enumerate all possible cones and to check their sta- 
bility at an overall cost of 0{N'^). Additional efforts 
can even bring the final complexity to 0{N'^ log(A^)) 
i.e. faster than the midpoint algorithm. This is il- 
lustrated on Fig. [2]where we observe that in practice 
SISCone runs faster than the typical implementations 
of the midpoint algorithm without a seed threshold 
and at least as fast as when a 1 GeV seed threshold 
is used. 

This has been implemented [SI El [S] in a C++ code 
named SISCone (Seedless Infrared Safe Cone) which 
is the first cone algorithm to satisfy the SNOWMASS 
requirements, that is to be at the same time IR and 
coUinear safe, and to be fast enough to be used in 
experimental analysis. 



■ CDF midpoint (s=0 GeV) 
CDF midpoint (s=1 GeV) 
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Figure 2: Clustering time for SIS- 
Cone compared to typical implemen- 
tations of the midpoint algorithm 
and the anti-kt algorithm [5] . 
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Figure 3: Jets found by the iterative cone for a 3-particle event (left) and for the same event 
with a coUinear splitting (right). 



Anti-fet as a replacement for the iterative cone. As for the midpoint algorithm, 
we start with an event with three hard particles (see Fig. EJa)). When clustered with the 
iterative cone, one stable cone containing all particles is found, resulting in a 1-jet event. 
If we now split the hardest particle into two coUinear particles — a process that also has 
an infinite probability in perturbative QCD — as shown on Fig. [3jb), clustering with the 
iterative cone gives a first jet made of particle 1 plus the two coUinear ones, then a second 
jet with particle 3. This example proves that the iterative cone algorithm is coUinear unsafe. 

Quite surprisingly, we can find a solution to that problem by coming back to the class of 
the recombination algorithms. The distance measures introduced earlier can be written as 



.{k',%kZ){Ayf, 



■A4) 



P,[GeV] 



with p = 1 for the kt algorithm and p = for the Cambridge/ Aachen algorithm. We can 
then consider a third case, the one for which p = —1 and call it the anti-kt algorithm [H]. 
Obviously, this algorithm is IR and coUinear safe. Furthermore, since its implementation [S] 
is similar to the one of the kt algorithm, its speed will be similar too, which certainly makes 
it usable for experimental purposes as seen on Fig. [2l 

To understand the link with the itera- 
tive cone algorithm, we note from the defini- 
tion of the anti-fct distance that pairs involv- 
ing a hard particle will be given small dis- 
tances. This means that soft particles are re- 
combined with hard ones before recombining 
among themselves, resulting in regular, soft- 
resilient, hard jets. This is exactly the hall- 
mark of the iterative cone and it is in that 
respect that the anti-fct can be seen as an IR 
and coUinear safe replacement. 

To illustrate this property, we show in Fig. 
m the jets resulting from the clustering of an 
event made with a few hard particles and a 
large number of very soft ones uniformly dis- 
tributed. It is clear that the hardest jets are 
perfectly circular and that, in general, the 
boundaries between the jets are regular. 




Figure 4: Illustration of the regularity of the 
jets obtained with the anti-fct algorithm. 
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Observable 


1st miss cones at 


Last meaningful order 


Inclusive jet cross section 


NNLO 


NLO 


W/Z/H + 1 jet cross section 


NNLO 


NLO 


3 jet cross section 


NLO 


LO (NLO in NLOJet) 


W/Z/H + 2 jet cross sect. 


NLO 


LO (NLO in MCFM) 


jet masses in 3 jets 


LO 


none (LO in NLOJet) 



Table 1: Perturbative level at which IR or coUinear unsafely arises for various processes. 



Physical impact and discussion. As we have seen, the seeded approach to stable 
cone search suffers from problems w.r.t. perturbative QCD expansion: the algorithms with 
split-merge are IR unsafe, while the iterative cone (with progressive removal) is collinear 
unsafe. We have introduced SISCone as a natural replacement of the cone algorithms with 
split-merge like midpoint, and the anti-Zci algorithm as a candidate to replace the iterative 
cone. These new algorithms are both IR and collinear safe. 

The question one might ask is to what extend these IR and collinear safety issue are 
important in real measurements. Since the unsafely arises when one has 3 particles in a 
common vicinity, it becomes important at the order af or aEwQ^s of the perturbative series. 

Table [T] summarises for different physical pro- 
cesses, the order at which seeded algorithms stop to 
be valid. The main message we can get from that ta- 
ble is thus that, if we do not want theoretical efforts 
in precise QCD computations to be done in vain, the 
resort of an IR and collinear safe algorithm like SIS- 
Cone and the anti-fcf is fundamental. To illustrate 
the argument more quantitatively. Fig. [5] shows the 
relative difference, expected to be present at the LO 
of perturbative QCD, between SISCone and midpoint 
for the mass of the 2"*^ hardest jet in 3-jet events. 
Differences reaching up to 40% are observed, proving 
that an IR and collinear safe algotithm is mandatory. 
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Figure 5: Relative difference between 
midpoint and SISCone for the mass 
of the 2"^^ hardest jet in 3-jet events. 
The 2"*^ and 3"* jets arc imposed to 
be distant by at most 2R. 
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